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I study fluid flow at the interface between elastic solids with randomly rough surfaces. I use 
the contact mechanics model of Persson to take into account the elastic interaction between the 
solid walls and the Bruggeman effective medium theory to account for the influence of the disorder 
on the fluid flow. I calculate the flow tensor which determines the pressure flow factor and, e.g., 
the leak-rate of static seals. I show how the perturbation treatment of Tripp can be extended 
to arbitrary order in the ratio between the root-mean-square roughness amplitude and the average 
interfacial surface separation. I introduce a matrix -D(C), determined by the surface roughness power 
spectrum, which can be used to describe the anisotropy of the surface at any magnification (. I 
present results for the asymmetry factor 7(1^) (generalized Peklenik number) for grinded steel and 
sandblasted PMMA surfaces. 



1. Introduction 

The influence of surface roughness on fluid flow at the 
interface between solids in stationary or sliding contact is 
a topic of great importance both in Nature and Technol- 
ogy. Technological applications includes leakage of seals, 
mixed lubrication, and removal of water from the tire- 
road footprint. In Nature fluid removal (squeeze-out) is 
important for adhesion and grip between the tree frog or 
Gecko adhesive toe pads and the countersurface during 
raining, and for cell adhesion. 

Almost all surfaces in Nature and most surfaces of 
interest in Tribology have roughness on many different 
length scales, sometimes extending from atomic distances 
(~ 1 nm) to the macroscopic size of the system which 
could be of order ^ 1 cm. Often the roughness is fractal- 
like so that when a small region is magnified (in general 
with different magnification in the parallel and orthog- 
onal directions) it "looks the same" as the unmagnified 
surface. 

Most objects produced in engineering have some par- 
ticular macroscopic shape characterized by a radius of 
curvature (which may vary over the surface of the solid) 
e.g., the radius i? of a cylinder in an engine. In this case 
the surface may appear perfectly smooth to the naked eye 
but at short enough length scale, in general much smaller 
than R, the surface will exhibit strong irregularities (sur- 
face roughness). The surface roughness power spectrum 
C(q) of such as surface will exhibit a roll-off wavelength 
Ao << i? (related to the roll-off wavevector qo — 2tt/Xo) 
and will appear smooth (except for the macroscopic cur- 
vature R) on length scales much longer than Aq. In this 
case, when studying the fluid flow between two macro- 
scopic solids, one may replace the microscopic equations 
of fluid dynamics with effective equations describing the 
average fluid flow on length scales much larger than Aq, 
and which can be used to study, e.g., the lubrication of 
the cylinder in an engine. This approach of eliminating 
or integrating out short length scale degrees of freedom to 
obtain effective equations of motion which describes the 
long distance (or slow) behavior is a very general and 
powerful concept often used in Physics. 



In the context of fluid flow at the interface between 
closely spaced solids with surface roughness, Patir and 
Cheng [1], have showed how the Navier- Stokes equa- 
tions of fluid dynamics can be reduced to effective equa- 
tions of motion involving locally averaged fluid pressure 
and flow velocities. In the effective equation occur so 
called flow factors, which are functions of the locally av- 
eraged interfacial surface separation u. They showed how 
the flow factors can be determined by solving numeri- 
cally the fluid flow in small rectangular units with lin- 
ear size of order (or larger than) the roll-off wavelength 
Aq introduced above. However, with the present speed 
(and memory) limitations of computers fully converged 
solutions using this approach can only take into account 
roughness over two or at most three decades in length 
scale. In addition, Patir and Cheng did not include the 
long-range elastic deformations of the solid walls in the 
analysis. Later studies have attempted to include elas- 
tic deformation using the contact mechanics model of 
Greenwood- Williamson (GW)0], but it is now known 
that this theory (and other asperity contact models 4]) 
does not correctly describe contact mechanics because of 
the neglect of the long range elastic coupling between 
the asperity contact regions [1, lij. In particular, the re- 
lation between the average interfacial separation u and 
the squeezing pressure p, which is very important for the 
fluid flow problem, is incorrectly described by the GW 
model [the GW model predict asymptotically (for large 
u) p ^ exp(— au^), while the exact resultjJQ for ran- 
domly rough surfaces is p ^ exp(— few), where a and b 
are constants determined by the nature of the surface 
roughness] . 

The paper by Patir and Cheng was followed by many 
other studies of how to eliminating or integrate out the 
surface roughness in fluid flow problems (see, e.g., the 
work by Sahlin et al. JlOJ ) . Most of these theories involves 
solving numerically for the fluid flow in rectangular inter- 
facial units and, just as in the Patir and Cheng approach, 
cannot include roughness on more than ~ 2 decades in 
length scale. In addition, in some of the studies the mea- 
sured roughness topography must be "processed" in a 



non-trivial way in order to obey periodic boundary con- 
ditions (which is necessary for the Fast Fourier Transform 
method used in some of these studies). 

Tripp[llj has presented an analytical derivation of the 
flow factors for the case where the separation between 
the surfaces is so large that no direct solid-solid contact 
occurs. He obtained the flow factors to first order in 
{h?) /u^, where {h?) is the ensemble average of the square 
of the roughness amplitude and u the average surface 
separation. This result is of great conceptual importance, 
but of minor practical importance as the influence of the 
surface roughness on the fluid flow becomes important 
only when direct solid-solid contact occur. 

Many surfaces of practical importance have roughness 
with isotropic statistical properties, e.g., sandblasted sur- 
faces or surfaces coated with particles typically bound 
by a resin to an otherwise flat surface, e.g., sandpa- 
per surfaces. However some surfaces of engineering in- 
terest have surface roughness with anisotropic statisti- 
cal properties, e.g., surfaces which have been polished 
or grinded in one direction. The theories of Patir and 
ChenP, [3] and of Tripp 11 1 can be applied also to sur- 
faces with anisotropic statistical properties. The surface 
anisotropy is usually characterized by a single number, 
the so called Peklenik number 7, which is the ratio be- 
tween the decay length of the height-height correlation 
function {h{x, y)h{0, 0)) along the x and y-directions, i.e., 
7 = ^x/^y where 0)/i(0, 0)) = (/i(0, 0)/i(0, 0))/2 

and (/i(0,Cy)/i(0,0)) = (/i(0,0)/i(0,0))/2. Here it has 
been assumed that the x-axis is oriented along one of 
the principal direction of the anisotropic surface rough- 
ness. However, the anisotropy properties of a surface may 
depend on the resolution (or magnification) which is not 
taken into account in this picture. 

In this paper I will present a new approach to calculate 
the fluid flow at the interface between two elastic solids 
with randomly rough surfaces. The present treatment is 
based on a recently developed theory for calculating the 
leak rate of stationary seals [l^- The theory use the con- 
tact mechanics theory of Persson[l^ [iJl in combination 
with the Bruggeman effective medium theory to calculate 
the fluid conductivity tensor. In this pap er we will gener- 
alize the treatment presented in Ref. [12| to surfaces with 
random roughness with anisotropic statistical properties. 
We also introduce a generalized Peklenik number 7(C) 
which depends on the magnification 7. Thus the theory 
takes into account that the anisotropy properties of the 
surface roughness may depend on the magnification un- 
der which the surface is observed. We present results for 
how 7(0 depends on ( for a grinded steel surface stud- 
ied using Atomic Force Microscopy and Scanning Tun- 
neling Microscopy, and for a sandblasted PMMA surface 
studied using an optical technique. As an illustration 
we calculate the pressure flow factor for surfaces with 
anisotropic properties. We emphasize that the present 
treatment accurately accounts for surface roughness on 
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FIG. 1: An elastic solid witli a rough surface in contact with 
a rigid solid with a flat surface. 



arbitrary many decades in length scale, and a full cal- 
culation typically takes less than a minute on a normal 
PC. In particular, the presented theory should be very 
useful for gaining a quick insight into what are the most 
important length scales in the problem under study. 

This paper is organized as follows: In Sec. 2 we briefly 
review the basic equations of fluid dynamics and describe 
some simplifications which are valid in the present case. 
In Sec. 3 and Appendix A I show how the perturbation 
treatment of Tripp can be extended to arbitrary order 
in {h'^)/u'^. This treatment may not be so important for 
the fluid fiow problem we consider as it is necessary to 
take into account that asperity contact occur already for 
relative small values of (ft.^)/u^, but the approach may 
find applications in other contexts. In addition, the solu- 
tion we present in wavevector space differ from the treat- 
ment of Tripp and leads directly to a matrix -D(C) which 
we used to describe the anisotropy of the surface at any 
magnification ^. In Sec. 4 we define I?(C) and present 
results for how the asymmetry factor 7(C) (generalized 
Peklenik number) depends on the magnification In 
Sec. 5 we briefly review the contact mechanics model we 
use. In Sec. 6 we describe the critical junction theory 
for the flow factor, and in Sec. 7 and 8 we show how the 
Bruggeman effective medium theory can be used in com- 
bination with the contact mechanics theory to calculate 
the fluid fiow tensor which determines the pressure flow 
factor and, e.g., the leak-rate of stationary seals. Sec. 9 
contains the summary. 

2. Fluid flow between solids with random sur- 
face roughness 

Consider two elastic solids with randomly rough sur- 
faces. Even if the solids are squeezed in contact, be- 
cause of the surface roughness there will in general be 
non-contact regions at the interface and, if the squeezing 
force is not too large, there will exist non-contact chan- 
nels from one side to the other side of the nominal con- 
tact region. We consider now fluid flow at the interface 
between the solids. We assume that the fluid is Newto- 
nian and that the fluid velocity fleld v(x, t) satisfies the 



Navier-Stokes equation: 
V • Vv 



(9v 

at 



1. 



where v = rj/p is the kinetic viscosity and p the mass den- 
sity. For simphcity we will also assume an incompressible 
fluid so that 

V • V = 

We assume that the non-linear term v • Vv can be ne- 
glected (which correspond to small inertia and small 
Reynolds number) , which is usually the case in fluid flow 
between narrowly spaced solid walls. For simplicity we 
assume the lower solid to be rigid with a flat surface, 
while the upper solid is elastic with a rough surface. In- 
troduce a coordinate system xyz with the xy-plane in the 
surface of the lower solid and the z-axis pointing towards 
the upper solid, see Fig. [TJ The upper solid moves with 
the velocity vq parallel to the lower solid. Let u{x, y, t) be 
the separation between the solid walls and assume that 
the slope |Vii| << 1. We also assume that u/L « 1, 
where L is the linear size of the nominal contact region. 
In this case one expect that the fluid velocity varies slowly 
with the coordinates x and y as compared to the varia- 
tion in the orthogonal direction z. Assuming a slow time 
dependence the Navier Stokes equations reduces to 

Here and in what follows v ~ (vx,Vy), x — {x,y) and 
V = {dx,dy) are two-dimensional vectors. Note that 
Vz ~ and that p(x) is independent of z to a good ap- 
proximation. The solution to the equations above can be 
written as 

1 z 
V = —z{z - u{x.))Wp + —r^vo 
Zl] u(x) 

so that V = on the solid wall z = and v = vq for 
z = u(x). Integrating over z (from z = to z = 



gives the fluid flow vector 



-^Vp+-u(x)vo 



(1) 



Mass conservation demand that 



5u(x, t) 
dt 



V • J = 



where the interfacial separation u{'x.,t) is the volume of 
fluid per unit area. In this last equation we have allowed 
for a slow time dependence of u(x, t) as would be the 
case, e.g., during fluid squeeze-out from the interfacial 
region between two solids. However, in this paper we 
will only focus on the case where u is time independent 



so that V • J = 0. This case is relevant for, e.g., fluid 
leakage in stationary seals. 

3. Perturbation treatment 

Here we show how one can obtain an effective flow 
equation by integrating out the short-wavelength rough- 
ness. We first re-derive the (first order) expansion re- 
sult of Tripp in wavevector space. After that we present 
the results of a Renormalization Group type of treat- 
ment (the derivation is presented in Appendix A). The 
treatment presented here does not take into account the 
elastic interaction between the solid walls and is there- 
fore strictly valid only for large enough average wall- wall 
separation. 

Let u(x) = u+h{x.) denote the local surface separation, 
where u = (u) is the average separation ((..) stands for 
ensemble averaging), and ^(x) is the contribution from 
the surface roughness with (h) = 0. In this section we as- 
sume h/u << 1 and perform a perturbation expansion in 
the small parameter h/u. Let us write the fiuid pressure 
as 

p = Po+Pi+P2 + ••• 

where po is the pressure to zero order in h (so that pq — 
(po)), pi to first order in h and so on. The fiuid fiow 
current is given by 

J = Vp H — uv 

Thus to second order in h we get 



3u^h„, , 3uh^ „ 1, ,, 

-^^^Po+Pi)"^^Po+2^^ + '^)^ (2) 

The ensemble average of this equation gives 



(J) =-Y2^V(po+Pi+P2) 



, 3S(/i2) 1 

where we have used that (h) = 0. Using that 

V-J = 
we get from (2) to zero order in h: 

V% = 0. 

The first order contribution gives 
3u^ 



(3) 



1 



V • V/i = (4) 



We define 



^^^^^ " (2^ / '^^^ ^1 We" 



and similar for /i(x). Substituting these results in (4) 
gives 

Y^9V(q)-Y^/i(q)(«q)-Vpo+2V-(*q)Mq) =0 (5) 
or 

Pi(q) = ^/i(q)(*q) • Vpo - 4^v • (iq)/i(q) (6) 
Next, note that 
(/i(q)/i(q')) = / {h{p.)h{pd))^'^-+^^-' 

= TTTTI I d''x(fx' (/i(x-x')/i(0))e''»-"+''i'-"' 



jq- (x-x')+i(q'+q) -x' 



(27r)^ 



^ / d^a;^^^' (/i(x-x')/i(0))e 



(2 



q-x+i(q'+q)-x' 



(fx (/i(x)/i(0))e''i-^5(q + q') 



= C(q)5(q + q') 
Using this equation and (6) gives 

(Wpi) = j <fqd?q' (iq')(/i(q)Pi(q'))e*^''+'''^-" 



Substituting this result in (3) gives 



(7) 



A and B can be written as ^ = u^cj)^ and B = ucj)^^ with 



where p = {po + Pi + P2) , and where the 2x2 matrices 
A and B can be written 
the flow factor matrices 

<^n = 1 



l((/.^)-3/AC(q)SJ 



1 + 



^h^) 



(1-3D), 



and 



3 /"^2 Nqq 

■ M 9 C'(q)- 



3(/i2) 



D. 



Here we have deflned the 2x2 matrix 



(8) 



(9) 



D 



f£q C(q)qq/g^ 



where g'o is the smallest surface roughness wavevector. 
For roughness with isotropic statistical properties, Dij = 
1/2 in which case (8) and (9) becomes 



2 m2 



(10) 



In deriving (7) we have used that to order h"^ one can 
replace terms like h'^Vpo with h'^Vp. 

In the derivation above we calculated the pressure and 
shear flow factors to first order in {h?)/u'^. In principle it 
is possible to extend the perturbation expansion to cal- 
culate higher order terms in This will result 
in higher order correlation functions, e.g., (/i 1/12/13/14) 
(where hi = /i(qi) and so on), but if the surface is ran- 
domly rough then these higher order correlation functions 
can be decomposed into a sum of products of pair corre- 
lation functions, e.g., 

(/11/12/13/14) = {hih2){hshi)+{hihi){h2hi)+{hihi){h2hs) 

Thus, all terms in the perturbation expansion will only 
involve the pair correlation function C (q) . We empathize 
that this is the case only for randomly rough surfaces 
where the phase of the different plane-wave components 
in the Fourier decomposition of /i(x) are uncorrelated. 
However, already the calculation of the second order term 
in the expansion of the flow factors in (/i^)/u^ becomes 
very cumbersome. In Appendix A we present a much 
simple and more powerful approach, which is in the spirit 
of the Renormalization Group (RG) procedure. Thus we 
eliminate or integrate out the surface roughness compo- 
nents in steps and obtain a set of RG flow equations de- 
scribing how the effective fluid equation evolves as more 
and more of the surface roughness components are elim- 
inated. 

Assume that after eliminating all the surface rough- 
ness components with wavevector |q| = q > Cqo the fluid 
current takes the form 

J = -^A{n, OVP + \b{u, C)v (11) 

where A and B are 2x2 matrices. In Appendix A we 
show that A{u,() and B{u,() satisfies 



OA 



^A"{u,0-A'{u,OMA'{u,0 



(12) 



dB 



-B"iu,0-A'iu,OMB'iu,0 



d_ 



(13) 



where A' = dA/du and so on, and where the 2x2 matrix 
M ^ A^^ is defined in Appendix A. Here {h?)c^ is the 
mean of the square of the roughness amplitude including 
only the roughness components with wavevector q > (.qq 
which can be written as 



d?q C(q) 



(14) 



?>C9o 



If we assume that D{Q) (defined in Appendix A and in 
Sec. 4) is independent of C, it is easy to solve these 
equations using perturbation theory to arbitrary order 
in the surface roughness amplitude h. As an example, 
for random roughness with isotropic statistical properties 
one obtain to second order in {h?) /v? (see Appendix A): 



B 



1 - 



3(^ 
2 u2 



3 (fe^)c 
2 u2 




(15) 



(16) 



The terms to linear order in (/i^) in these expressions 
agree with the result of Tripp. He compared his ex- 
pansion results with the numerical results of Patir and 
Cheng and found that the expression for A (or 0p) and 
B (or (ps) agree rather well with the numerical results 
for (Jx^)^!"^ ju < 3 and < 6, respectively. For the lat- 
ter case our second order contribution to B improves the 
agreement between numerical results and the expansion 
result but for {h?)^/"^ /u < 3 the direct wall- wall inter- 
action becomes so important that the expansion result 
(which neglect this interaction) cannot be used. 

4. Surfaces with anisotropic statistical proper- 
ties 

As discussed in the introduction, surfaces with 
anisotropic statistical properties are usually character- 
ized by the Peklenik number 7 — ^x/^y, which is the 
ratio between characteristic correlation length and , 
defined as the distances along the x and y-axis where 
the height-height correlation function {h{x, y)h{0, 0)) has 
decayed to half of its initial value. However, for most 
real surfaces 7(C) will depend on the magnification or 
length-scale under consideration. Here we propose to 
obtain 7(C) from the surface roughness power spectrum 
C{qx,qy) as follows: 

The surface roughness power spectrum C(q) is defined 

by 



1 



(fx (/i(x)/i(0))e 



iq-x 



(a) AFM 




(b) STM 




FIG. 2: Surface topography of a grinded steel surface obtained 
using (a) Atomic Force Microscopy (AFM) (10/im x lO^m) 
and (b) Scanning Tunneling Microscopy (STM) (0.1/im x 
O.l^im). 
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FIG. 3: The (angular averaged) surface roughness power spec- 
trum C{q) calculated from the AFM and the STM surface 
topography data shown in Fig. [2] 



We can write 

C(x) = (/i(x)/i(0)) = j fq C(q)e-''i-" 
We also define 

2tt 



C(x,C) = / d0C(q) 
Jo 



where q = C(7o(cos0, sin0). Now consider the closed con- 



tour defined by 



C(x,C) =C(0,C)/2 



We now fit this contour to tlie quadratic function /(x) — 
+biXi+c. The function = const, describes 

an ellipse which in general has its major axis rotated 
by some angle ip relative to the x-axis. We define 7 as 
the ratio between the major and minor ellipse axis, and 
obtain both 7(^) and the rotation angle V'(C)i both of 
which depend on the magnification 

Another way to determine an effective 7(C) is as fol- 
lows: Consider the tensor (see also Appendix A) 



^(0 



C(q)qq/q2 



/o 



C(q) 



(17a) 



where q = ^90 (cos0, sin(/)) . If D{() is independent of C 
then this definition is identical to 



D 



Jd^q C(q)qq/g^ 
Jd^q C(q) 



(17b) 



which appeared already in the perturbation calculation 
in Sec. 3. Note that Dn + D22 = TrD ^ 1 and that the 
D is symmetric and can be diagonalized. For example, 
suppose C(q) — f{axql + ciyqy) and that the q-integrals 
in (17b) are over the whole q-plane. For this case we get 
after some simplifications 



D 



1 

2^ 



where 7^ = ax/oty- Performing the integral gives Dn = 
1/(1 +7) and D22 = 7/(1 + 7). Note that in this case 
\D\ = D11D22 = 7/(1 + 7)^ where \D\ is the determinant 
of the matrix D. This equation has two solutions, 7 and 
1/7 where 



1 



7 



2\D\ 



i-n-m) 



1/2 



- 1 



(19) 



Note that this definition of 7 is independent of the coordi- 
nate system used since the determinant is invariant under 
rotations (orthogonal transformations). Note also that 
for a surface with isotropic statistical properties from 
(17) = S,j/2 so that \D\ = 1/4 and (19) reduces 
to 7 = 1 as it should. The angle ?/> between the major 
axis of the ellipse and the x-axis of the coordinate sys- 
tem depends, of course, on the coordinate system, and is 
given by 



tarn/) = c ± (1 



,1/2 



(20) 



where c = {D22 - Dii)/{2Di2)- 

In Fig. [2] we show the surface topography of a grinded 
steel surface as obtained using (a) Atomic Force Mi- 
croscopy (AFM) (10/zm X lO^m) and (b) Scanning Tun- 
neling Microscopy (STM) (0.1/xm x 0.1/^m). In Fig. [3] 
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FIG. 4: The 7-parameter calculated from the AFM and the 
STM surface topography data shown in Fig. [2] The maxi- 
mum of 7 occur for q ~ 1.8 x 10® m~^, corresponding to a 
wavelength A = 2n/q ~ 3.5 ^m. This is just the wavelength 
of surface topography orthogonal to the major wear tracks in 
Fig. 



I show the (angular averaged) surface roughness power 
spectrum C{q) calculated from the AFM and the STM 
surface topography data shown in Fig. [2] The power 
spectrum is well approximated with self affine fractal 
with the fractal dimension Df — 2.25. However, note that 
the surface topography is anisotropic. In Fig. |4]we show 
the calculated (using (19)) 7-parameter for the same sur- 



(18) face. The maximum of 7 occur for q w 1.8 x 10^ 



corresponding to a wavelength A = 2'!T/q w 3.5 /xm. This 
is just the wavelength of the surface topography orthog- 
onal to the major wear tracks in Fig. [2] 

In Fig. [5] we show the calculated (using (19)) 7- 
parameter for a sandblasted PMMA surface. In this case 
the statistical properties of the surface are expected to 
be isotropic, and indeed 7 is very close to unity. 

For surfaces which have been grinded or polished in one 
direction, wear scars may occur almost uninterrupted for 
a very long distance. In this case it is necessary to mea- 
sure the surface topography over a very large surface area 
in order to correctly obtain the 7(C)-function. In numer- 
ical fiow calculations as involved in, e.g., the studies of 
Patir and Cheng, it would be necessary to use very large 
rectangular units which would be practically impossible 
because of the huge memory and computational time re- 
quired. 

5. Contact mechanics: short review and basic 
equations 

At short (average) interfacial separation there will be a 
direct asperity interaction between the solids walls, and 
in this case the perturbation approach of Sec. 2 will fail. 
Here we will briefly review the contact mechanics model 
of Persson which we use in this study. 
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FIG. 5: The 7-parameter calculated from optically measured 
surface topography data for sandblasted PMMA. The surface 
topography was measured over a 3 cm x 3 cm surface area. 
The surface root-mean-squarc roughness was 32 jum. 



magnification C, 



elastic sotid 




FIG. 7: An asperity contact region observed at the magnifica- 
tion (,. It appears that complete contact occur in the asperity 
contact region, but when the magnification is increasing to the 
highest (atomic scale) magnification ^1 , it is observed that the 
solids are actually separated by the average distance w(C). 




FIG. 6: An rubber block (dotted area) in adhesive contact 
with a hard rough substrate (dashed area). The substrate 
has roughness on many different length scales and the rubber 
makes partial contact with the substrate on all length scales. 
When a contact area is studied at low magnification it appears 
as if complete contact occur, but when the magnification is 
increased it is observed that in reality only partial contact 
occur. 



Consider the frictionless contact between two elastic 
solids with the Young's elastic modulus Eq and Ei and 
the Poisson ratios and vi . Assume that the soUd sur- 
faces have the height profiles /io(x) and /ii(x), respec- 
tively. The elastic contact mechanics for the solids is 
equivalent to those of a rigid substrate with the height 
profile /i(x) = ho{x) + /ii(x) and a second elastic sohd 
with a flat surface and with the Young's modulus E and 
the Poisson ratio i' chosen so that 15 1 



The contact mechanics formalism developed 
elsewhere 0, i, [H, Eil is based on the studying the 
interface between two contacting solids at different 
magnification C. When the system is studied at the mag- 
nification C it appears as if the contact area (projected on 
the xy-plane) equals A(C), but when the magnification 
increases it is observed that the contact is incomplete 
(see Fig. and the surfaces in the apparent contact 
area A{() are in fact separated by the average distance 
u{C), see Fig. [T) The (apparent) relative contact area 
A{(^)/Aq at the magnification C is given bvjsl. [l3| 

Ao " (vrG)V2 ' 
where 



erf I 



G(C) 



E 



dqq'Ciq) 



where the surface roughness power spectrum 



Ciq) 



(fx (/i(x)/i(0))e" 



iq-x 



where (...) stands for ensemble average. The height pro- 
file /i(x) of the rough surface can be measured routinely 
today on all relevant length scales using optical and sty- 
lus experiments. 

The quantity u{C.) is the average separation between 
the surfaces in the apparent contact regions observed 
at the magnification see Fig. [71 It can be calculated 
fromfsj 

"(0 = %/^ Trf? q'C{q)w{qX) 
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FIG. 8: (a) The black area is the asperity contact regions at 
the magnification The green area is the additional contact 
area observed when the magnification is reduced to ( — A( 
(where Ai^ is small) . The average separation between the solid 
walls in the green surface area is denoted by (b) The 

separation between the solid walls along the blue dashed line 
in (a). Since the surfaces of the solids are everywhere rough 
the actual separation between the solid walls in the green area 
will fluctuate around the average At the most narrow 

constriction the surface separation is u^. 



where — poAo/A{Q and 



w{q,C) 



q X -1/2 
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We define ui(C) to be the (average) height separating 
the surfaces which appear to come into contact when the 
magnification decreases from ^ to C — A^, where AC is a 
small (infinitesimal) change in the magnification. In Fig. 
[5]Ja) the black area is the asperity contact regions at the 
magnification C- The green area is the additional contact 
area observed when the magnification is reduced to C~ A^ 
(where AC is small) [l^ . The average separation between 
the solid walls in the green surface area is given by ui(C)- 
Fig. |Sl^b) shows the separation between the solid walls 
along the dashed line in Fig. Since the surfaces 

of the solids are everywhere rough the actual separation 
between the solid walls in the green area will fluctuate 
around the average wi(C)- Thus we expect the smallest 
surface separation Uc = aui(Cc)^ where a < 1 (but of 
order unity, see Fig. [8i;b))0. In Ref. [H we have 
analyzed leak-rate data for rubber seals and always found 
that a to be in the range 0.5 — 1. However, it is clear 
that a cannot be a fixed constant but must depend on the 



average surface separation and on the surface roughness 
which occur at length scales shorter than A = L/C. In 
particular, as {h^) / u\{C} — >■ we expect that a — > 1 (see 
also Sec. 8). 

ui(C) is a monotonically decreasing function of C, and 
can be calculated from the average interfacial separation 
m(C) and A(C) using (see Ref. ^) 

«i(C)=«(C) + ^'(C)A(C)M'(C)- 

One can show[l3] from the equations above that as the 
applied squeezing pressure po — > 0, for the magnifications 
most relevant for calculating fluid flow (e.g., the leak-rate 
of seals). Ml — u. 

We note that when solving for the fluid flow between 
macroscopic surfaces with roughness one may in a mean- 
field type of treatment write the local nominal pressure 
(i.e., the pressure locally averaged over surface area with 
linear dimension of order the wavelength Aq of the longest 
surface roughness component) as(l9| 

t) = Pfluid(x, t) + Psolid(x, t) 

where pfiuid and Psoiid are locally averaged nominal fluid 
pressure and solid wall-wall contact pressure, respec- 
tively. The pressure Psoiid can be related to the inter- 
facial separation m(x, t) as described in Ref. 0, HI- In 
particular, for large enough average surface separation^ 

Psoiid ~ /3£;*e-"/"° 

where /? and can be calculated from the surface rough- 
ness power spectrum. 

6. Critical-junction theory of fluid flovif 

The perturbation expansion presented in Sec. 3 as- 
sumed no direct contact between the solid walls. But di- 
rect contact between the solid walls occur in most cases 
of interest, e.g., in static seals. The simplest approach 
for this case is based on the leak-rate model developed 
m Ref. [i2,[23-l23. Consider the fluid leakage through a 
(nominal) contact region, say between a hard solid and 
rubber, from a high fluid pressure Pa region, to a low fluid 
pressure Pb region. Assume that the nominal contact re- 
gion between the rubber and the hard countersurface is 
rectangular with area L^xLy, with Ly > L^. We assume 
that the high pressure fluid region is for a; < and the 
low pressure region for x > L^- We "divide" the contact 
region into squares with the side — L and the area 
Aq = (this assumes that N = Ly/Lx is an integer, but 
this restriction does not affect the final result). Now, let 
us study the contact between the two solids within one of 
the squares as we change the magnification C. We define 
C — i/A, where A is the resolution. We study how the 
apparent contact area (projected on the xy-plane), A(C), 
between the two solids depends on the magnification C- 
At the lowest magnification we cannot observe any sur- 
face roughness, and the contact between the solids ap- 
pears to be complete i.e., A{1) = Aq. As we increase the 
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critical constriction 



FIG. 9: The contact region at different magnifications C = 3, 
9, 12 and 648, is shown in (a)-(d) respectively. When the 
magnification increases from 9 to 12 the non-contact region 
percolate. At the lowest magnification!^ — 1: A{1) = Aq. The 
figure is the result of Molecular Dynamics simulations of the 
contact between elastic solids with randomly rough surfaces, 
see Ref. [H. 



magnification we will observe some interfacial roughness, 
and the (apparent) contact area will decrease. At high 
enough magnification, say C = Cci a percolating path of 
non-contact area will be observed for the first time, see 
Fig. ini We denote the most narrow constriction along 
this percolation path as the critical constriction. The 
critical constriction will have the lateral size Ac — i/Cc 
and the surface separation at this point is denoted by 
Uc = aui(^c)- As we continue to increase the magnifica- 
tion we will find more percolating channels between the 
surfaces, but these will have more narrow constrictions 
than the first channel which appears at C = Cc, and as a 
first approximation one may neglect the contribution to 
the leak-rate from these channels [2^. 

A first rough estimate of the leak-rate is obtained by 
assuming that all the leakage occurs through the critical 
percolation channel, and that the whole pressure drop 
AP = Pa — ft (where Pa and Pb is the pressure to the 
left and right of the seal) occurs over the critical con- 
striction (of width and length Ac ~ i/Cc and height Mc)- 
We will refer to this theory as the "critical-junction" the- 
ory. If we approximate the critical constriction as a pore 
with rectangular cross section (width and length Ac and 
height Uc << Ac), and if we assume an incompressible 
Newtonian fluid, the volume-fiow per unit time through 



the critical constriction will be given by (Poiseuille flow) 
Q ^ I^AP, (22) 

where rj is the fluid viscosity. In deriving (22) we have 
assumed laminar flow and that Uf. << Ac, which is al- 
ways satisfied in practice. We have also assumed no-slip 
boundary condition on the solid walls. This assumption 
is not always satisfied at the micro or nano-scale, but is 
likely to be a very good approximation in the present case 
owing to surface roughness which occurs at length-scales 
shorter than the size of the critical constriction. Finally, 
since there are N = Ly/Lx square areas in the rubber- 
countersurface (apparent) contact area, we get the total 
leak-rate 

,3 

-AP. (23) 



^ 12ri' 



Note that a given percolation channel could have several 
narrow (critical or nearly critical) constrictions of nearly 
the same dimension which would reduce the fiow along 
the channel. But in this case one would also expect more 
channels from the high to the low fluid pressure side of 
the junction, which would tend to increase the leak rate. 
These two effects will, at least in the simplest picture 
where one assumes that the distance between the critical 
junctions along a percolation path (in the x-direction) is 
the same as the distance between the percolation chan- 
nels (in the y-direction), compensate each other (see Ref. 
(2^). The effective medium theory presented below in- 
cludes (in an approximate way) all the flow channels. 

To complete the theory we must calculate the separa- 
tion Uc of the surfaces at the critical constriction. We 
first determine the critical magnification ^c by assum- 
ing that the apparent relative contact area at this point 
is given by percolation theory. Thus, the relative con- 
tact area A{C)/Ao ss 1 — pc, where Pc is the so called 
percolation threshold [lH. For infinite-sized 2D systems, 
and assuming site percolation, pc ~ 0.70 for a hexagonal 
lattice, 0.59 for a square lattice, and 0.5 for a triangu- 
lar latticefi^. For bond percolation the corresponding 
numbers are 0.65, 0.5, and 0.35, respectively. For contin- 
uous percolation in 2D the Bruggeman effective medium 
theory predict pc — 0.5. For finite sized systems the per- 
colation will, on the average, occur for (slightly) smaller 
values of Pc, and fluctuations in the percolation thresh- 
old will occur between different realizations of the same 
physical system. Numerical simulations such as those 
presented in Ref. [13] (see Fig. [9|) and Ref. [iij typically 
gives Pc slightly larger than 0.5. In our earlier leak-rate 
studies we have used Pc = 0.5 and 0.6 to determine the 
critical magnification — (^c- 

We can write the leak-rate in terms of the pressure flow 
factor. Thus the current 
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FIG. 10: Effective medium theories take into account random 
disorder in a piiysical system, e.g., random fluctuations in the 
interfacial separation m(x). Thus, for a n-component system 
(e.g., where the separation u takes n different discrete values) 
the flow in the effective medium should be the same as the 
average fluid flow obtained when circular regions of the n- 
components are embedded in the effective medium. Thus, 
for example, the pressure p at the origin calculated assuming 
that the effective medium occur everywhere must equal the 
average ^ dpi (where Ci is the concentration of component i) 
of the pressures pi (at the origin) calculated with the circular 
inclusion of component i = 1, n. 



and the leak-rate 



12?7 

Comparing this with (23) gives 

wi(C< 



a- 



u{l) 



7. Effective medium theory of fluid flow: 
isotropic roughness 

The critical-junction theory presented above assumes 
that the leak-rate is determined by the resistance towards 
fluid flow through the critical constriction (or through 
a network of critical constrictions, see above). In real- 
ity there will be many flow channels at the interface. 
Here we will use the 2D Bruggeman effective medium 
theory [23, [2^ to calculate (approximately) the leak-rate 
resulting from the network of flow channels. Another ap- 
proach to extend the critical junction theory is critical 
path analysis, see Ref. (29l.[30|. 

We study the fluid flow through an interface where 
the separation m(x) between the surfaces varies with the 
lateral coordinate x = {x,y). If u(x) varies slowly with 
X the Navier-Stokes equations of fluid flow reduces to 



J = -crVp 



(24) 



where the conductivity a — u^{-x)/12ri. 

In the effective medium approach one replace the lo- 
cal, spatial varying, conductivity cr(x) with a constant 
effective conductivity a^s- Thus the fluid flow current 
equation 



as applied to a rectangular region x Ly with the pres- 
sure gradient dp/dx = (Pb — Pa)/Lx, gives 



Ly 

LyJx =^ — CTcff AP 



(26) 



where AP — Pn — Ph is the pressure drop. 

The effective medium conductivity doff is obtained as 
follows. Let us study the current flow at a circular in- 
clusion (radius R) with the (constant) conductivity a lo- 
cated in an infinite conducting sheet with the (constant) 
conductivity ficff- We introduce polar coordinates with 
the origin at the center of the circular inclusion. The 
current 

J = — trVp for r < R 

J = — (JoffVp for r > R 
We consider a steady state so that 
VJ = 



or 



V^p = 



(27) 



If Jo = —(Toff a is the current far from the inclusion (as- 
sumed to be constant) we get for r > R: 

p= [l + /(r)]a.x (28) 

Eq. (27) is satisfied if 

/"(r) + 3/'(r)r-i =0 

A solution to this equation is / = ar^^. Substituting 
this in (28) gives 

p= [l + ar^^] a-x (29) 

For r < R we have the solution 

p = /3a • X (30) 

Since p and x • J must be continuous at r = i? we get 
from (28) and (29): 

1 -I- ai?~2 ^ p 

(1 — ai?^^) (Tcff — fHa 
Combining these two equations gives 

r, 2o-off 



CTcff + cr 



(31) 



(25) 



The basic picture behind effective medium theories is pre- 
sented in Fig. 1101 Thus, for a two component system, 
one assumes that the flow in the effective medium should 



be the same as the average fluid flow obtained when cir- 
cular regions of the two components are embedded in 
the effective medium. Thus, for example, the pressure 
p calculated assuming that the effective medium occur 
everywhere must equal the average cipi + C2P2 of the 
pressures pi and p2 calculated with the circular inclusion 
of the two components 1 and 2, respectively. For r < R 
we have for the effective medium p = a - x and using (30) 
the equation p = c\p\ + C2P2 gives 



1 = ci/3i + C2P2 



(32) 



where ci and C2 are the fractions of the total area oc- 
cupied by the components 1 and 2, respectively. Using 
(31) and (32) gives 



1 



ci- 



2q'eff 

fJeff + (Ti 



C2 



2(7eff 
fJeff + (72 



which is the standard Bruggcman effective medium for a 
two component system. Note that if one component is 
insulating, say (J2 = 0, as ci — )• 0.5 from above, Ueff 0, 
i.e., = 1/2 is the percolation threshold of the two 
component 2D-Bruggeman effective medium model. 

If one instead have a continuous distribution of com- 
ponents (which we number by the continuous index ^) 
with conductivities a = (t(^), then 



= jd^ p{om 



(33) 



where P{^) is the fraction of the total surface area oc- 
cupied by the component denoted by ^. The probability 
distribution P(^) is normalized so that 



/ 



p{0 = 1 



Using (31) we get 



2(Teff 



(7eff + O-(0 



(34) 



(35) 



It is easy to show from this equation that also for the 
case of a continuous distribution of components, the per- 
colation limit occur when the non-conducting component 
(which in our case correspond to the area of real contact 
where u = and hence a = /\2ri = 0) occupies 50% of 
the total surface area, i.e., Pc = 1/2 in this case too. 

To summarize, using the 2D Bruggeman effective 
medium theory we get: 



(a) 



(b) 



(0) 



5x 



•i I 



FIG. 11: Contact regions for (a) longitudinal oriented, (b) 
isotropic, and (c) transversely oriented rough surfaces. The 
ratio between the ellipse major a^xis is denoted by 7 = ^x/^y 
and 7 > 1, = 1 and < 1 in (a), (b) and (c), respectively. The 
average fluid flow is in the a;-direction. 
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FIG. 12: The black area denote interfacial solid-solid contact 
with the flow conductivity a2 = 0. The two cases (a) and (b) 
correspond to 7 = 00 and 7 = 0, respectively. In the first 
case (a) fluid flow can occur in the strips (open channels) of 
component 1 for arbitrary low concentration of component 1. 
In this case fluid flow wiU occur at the interface until complete 
contact occur between the solids. In the opposite limit 7 — >■ 
no fluid can flow (in the a-direction) at the interface unless 
C2 is zero. 
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where 



^(0 



Ao J aes + o-(C) ' 
12r] ■ 



(37) 



(38) 



Eq. (37) is easy to solve by iteration. 

8. Effective medium theory of fluid flow: 

anisotropic roughness 

Here we briefly describe how one may apply the effec- 
tive medium theory to study fluid flow between surfaces 
with anisotropic (but translational invariant) statistical 
properties. Let p be the locally averages pressure and J 
the fluid flow current also locally averaged. We have 



,e dp 
'■' dxi 



(39) 



(36) 



where AP = — -Pb is the pressure drop and where 

2 



— =[ da P{a) 



feff + cr 



Note that 



,eff 



We can choose a coordinate system such that the flow 
conductivity tensor is diagonal: 



CTcff = 



'^11 





In this case the x and y-coordinate axis are oriented 
along and perpendicular to the "groves" on the surface, 
respectively. The flow conductivity for any other orienta- 
tion can be obtained using the standard transformation 
of tensors under rotation. Thus if the x axis is oriented 
an angle relative to the "groves" then 



Ceff 



(T||cos^0 + crj_sin^0 (cr|| — (T^)cos(/)sin(/ 
((T|| — o'j_)cos0sin(/< (T||Sin^(/) + (T^cos^t; 



We will now calculate the flow conductivities and 
(T_L parallel and perpendicular to the groves, respectively. 
We assume that interfacial separation u(x) varies slowly 
with X — (x, y). Consider an elliptic inclusion in a fluid. 
Assume that the fluid flow conductivity equals cTcff out- 
side the inclusion and ui in the inclusion. Assume that 
the fluid flow far from the inclusion is oriented at an an- 
gle (f) relative to a major axis of the inclusion, i.e., far 
from the inclusion J = — fToff a and 



P 



a{xcos(j) + ysinq 



(40) 



The fluid flow can be calculated analytically using ellip- 
tic coordinates (^, i?), see Ref. [31j]. In this coordinate 
system the curves /i = const, are ellipses. Consider the 
ellipse fi = fj,o. The ratio 7 between the major and minor 
axis can be written as 7 = coth/iQ so that when 00 
the ellipse becomes a circle. 

The fluid pressure inside the elliptic inclusion is given 

by 



(41) 



where the matrix Aij has the components A12 — A21 — 
and 



An 



o'cfie'' 



22 



CTcff cosh/io + CTisinh/io 



o-cffe'^" 
CTiCosh/io -f- CTcffSinh/iQ 



(42) 



(43) 



Note that when /iq — ^ 00 the matrix Aij = pSij, where 
/3 is given by (A8). Thus in this case the pressure in the 
inclusion becomes p — /?a-x just as for a circular inclusion 
(see Eq. (A7)), which of course is expected because the 
ellipse becomes a circle when ^0 — ^ 00. 

The parameter ^0 is determined by 7 = coth/iQ or 



7 + 1 
7-1 



Using this equation we can also write (42) and (43) as 

o-off(7 + 1) 



An 



CTcff7 + CTi 



(44) 



^22 = 



r^csil + 1) 



(45) 



0-17-1- o-off 

We now consider the situation where (p = so that one 
of the ellipse axis is oriented along the (average) fluid flow 
direction as in Fig. [TTJa). In this case the pressure in 
the inclusion 

p — Anax 

while the pressure far away from the inclusion p = ax. 
For a two-component system the effective medium equa- 
tion (33) now becomes 



ClCTeff(7l + 1) C2CTcff (72 + 1) 



f off 71 + ^1 



0'off72 + f-2 



(46) 



where we have taken into account that the two compo- 
nents may have different ratio 7. Assume that one com- 
ponent, say component 2, has the conductivity CT2 = 0. 
In this case it follows from (46) that CTcff — ?> as C2 — > 
71/(71 + 1)- Note in particular that for 71 — 7> 00, C2 — > 1 
i.e. in the limit when the major axis of the inclusion 1 
goes to infinite (where conducting strips of the conduct- 
ing component 1 occur for arbitrary low concentration of 
component 1) fluid flow will occur at the interface until 
complete contact occur between the solids. In the oppo- 
site limit 71 — 0, C2 — >■ 0. In this case no fluid can flow 
(in the x-direction) at the interface for any applied pres- 
sure. These two limits correspond to the configurations 
illustrated in Fig. [TH 

For a continuous distribution of components 



1 - / rfC PiO 



gcff(7(C) + 1) 
f^cff7(C) + ct(C) 



(47) 



where ct(C) = (aMi(C))^/12?7. This equation is also valid 
for the orientation of the ellipse as in Fig. [TlT c) in which 
case 7 < 1 (in general, 7 is the ratio between the ellipse 
axis in the x-direction and the y-direction) . 

In Fig. [13] we show the pressure flow factor (f>p as a 
function of the average surface separation u in units of 
the root-mean-square roughness amplitude. In the cal- 
culation we have for simplicity assumed that 7(C) is a 
constant independent of the magnification C. Results are 
shown for three different surfaces with surface roughness 
with isotropic statistical properties (7 = 1), and for sur- 
faces with anisotropic roughness of longitudinal (7 = 3) 
and transverse (7 = 1/3) type. The dashed lines is calcu- 
lated with a = 1 while the solid lines is calculated with a 
a which depends on the interfacial separation as follows: 

As pointed out is Sec. 5, the surfaces in the (non- 
contact) flow channels are everywhere rough, and the ac- 
tual separation between the solid walls in the non-contact 
region which appears when the magnification is reduced 
from C to C — AC (green area in Fig. (H^a)) will fiuc- 
tuate around the average wi(C)- Thus with respect to 
fluid flow the separation u between the walls will ap- 
pear smaller than the average ui and we use u = aui , 




2 3 4 

average separation / rms 

FIG. 13: The pressure flow factor <j)p as a function of the 
average surface separation u in units of the root-mean-square 
roughness amplitude. For three different surfaces with surface 
roughness with isotropic statistical properties (7 = 1), and 
for surfaces with anisotropic roughness of longitudinal (7 = 
3) and transverse (7 = 1/3) type. The 7=1 case is for 
sandblasted PMMA (root-mean-square roughness 22 ^m) in 
contact with rubber with the elastic modulus E — 2.3 MPa. 
The other cases assumes the same angular averaged power 
spectrum and elastic properties as for the 7 = 1 case. The 
solid and dashed lines are discussed in the text. 
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FIG. 14: The variation of the area of real contact A (in units 
of the nominal contact area Aq) and the average interfacial 
separation u (in units of the root-mean-square roughness am- 
plitude) as a function of the (nominal) squeezing pressure 
for the system shown in Fig. 1131 sandblasted PMMA (root- 
mean-square roughness 22 /im) in contact with rubber with 
the elastic modulus E = 2.3 MPa. 



where a < 1 . We note that a is due to the surface rough- 
ness which occur at length scales shorter than A = i/C, 
and it is possible to calculate (or estimate) a from the 
surface roughness power spectrum, as follows. 

As shown in Sec. 2 and Appendix A, the fluid flow be- 
tween two nominal flat surfaces is affected by the surface 
roughness on the solid walls even at such large (average) 
surface separation that no direct wall- wall contact occur. 



Thus for isotropic roughness at large separation there is 
a reduction in the fluid flow entering via the flow factor 
(f)p ^ I ~ (3/2)((/i^)/m^), where u is the average surface 
separation. If we apply this to the present case in the 
fluid flow problem we replace the term (C) by (00* 
where 



3(/^2) 



Here we have assumed surface roughness with isotropic 
statistical properties and {h'^)(; denote the ensemble av- 
erage of the square of the roughness amplitude including 
only the surface roughness with wavevectors larger than 
q = C*'?o- In calculating the solid lines in Fig. [T3]we have 
chosen (* = 3C. 

Fig. [T51 shows, as expected, that when 7 decreases the 
percolation limit, below which no fluid flow can occur, ap- 
pears at larger and larger average separation. Note also 
that for 7 = 3 the pressure flow factor flrst increases with 
decreasing u, but finally it decreases towards zero. Thus, 
even for arbitrary large 7 at high enough squeezing pres- 
sures (corresponding to small enough u) the non-contact 
area will not percolate in which case no fluid flow can 
occur at the interface and (f)p — 0. 

Fig. [13] shows the variation of the area of real con- 
tact A (in units of the nominal contact area ^0) and 
the average interfacial separation u (in units of the root- 
mean-square roughness amplitude) as a function of the 
(nominal) squeezing pressure for the system shown in 
Fig. [131 sandblasted PMMA (root-mean-square rough- 
ness 22 fim) in contact with rubber with the elastic modu- 
lus E — 2.3 MPa. Note that even at the lowest squeezing 
pressure where u/rms « 4 the area of real contact is still 
non-negligible, about 1% of the nominal contact area. 

9. Summary and conclusion 

I have studied the fluid flow at the interface between 
elastic solids with randomly rough surfaces. I have used 
the contact mechanics model of Persson to take into ac- 
count the elastic interaction between the solid walls and 
the Bruggeman effective medium theory to account for 
the influence of the disorder on the fluid flow. I have cal- 
culate the flow tensor which determines the pressure flow 
factor and, e.g., the leak-rate of seals. I have shown how 
the perturbation treatment of Tripp can be extended to 
arbitrary order in the ratio between the root-mean-square 
roughness amplitude and the average interfacial surface 
separation. I have introduced a matrix I?(C), determined 
by the surface roughness power spectrum, which can be 
used to describe the anisotropy of the surface at any mag- 
niflcation Q. I have present results for the asymmetry 
factor 7(0 (generalized Peklenik number) for a grinded 
steel surface and a sandblasted PMMA surface. 
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Appendix A 

In Sec. 2 wc calciilatcd the prcssiirc and shear flow fac- 
tors to first order in {h^)/u^. Here we will present a sim- 
pler and more powerful approach, which is in the spirit 
of the Rcnormalization Group (RG) procedure. Thus 
we will eliminate or integrate out the surface roughness 
components in steps and obtain a set of RG flow equa- 
tions describing how the effective fluid equation evolves 
as more and more of the surface roughness components 
are eliminated. 

Assume that after eliminating all the surface rough- 
ness components with wavevector |q| = g > C5o the fluid 
current [given by (1)] takes the form 

3 = -±-A{u)Vp+\B{u)y; {Al) 

where A and B are 2x2 matrices. We now add to u a 
small amount of roughness 



h= dPqh{ci)e''^-'' 

•^«-AC)go<</<C9o 

Consider now the current 

J = - t^^(m + h)Wp + \b{u + h)v 
Writing as before 

we get to second order in h 



j = 



12r; 



■V(po +Pi +P2) 



-V(po hr—^VQ 



12r] 



2477 



+ + B'{u)h)v + ^B"{u)h^v {A3) 

The ensemble average of this current gives to second or- 
der in h 



^ -{hVpi) ' Vp 



+^B{u)^r+-^B"{u){h')v 



where we have used that (h) = 0. To zero order in h the 
continuity equation V • J gives 

Aij{u)didjPo = 0, 

and to first order in h we get 

_^d,hdjpo - ^d,djp, + \B',,{u)d,hvj = 

In wavevector space this equation takes the form 

1 „ J. 

12r] 



-r^Aj{u){iqi)h{q)djPo + —Aij{u)qiqjPi{ci) 



+-Bij{u){iqi)h{ci)vj=0 



or 



1277 



2477 



Pi{q) = {Aim{u)qiqm) ^{iqi)h{q) 

x{A',^{u)djpo-6rjBljiu)vj) {A5) 
Using this equation and (A2) gives 

{hdip,} = j d\d\'{h{c{){iqi)pM)) 

= J (fq C{q){Aim{u)qiqm)~^qiqj 

X {6r]B'.^{u)vk - A'.^{u)dkPo) {A6) 
Let us define the matrix 

My = j <fq C{q){Aim{u)qiqm)-\iqj {A7) 

so that (A5) becomes 

(hVpi) = 6r]MB'v - MA'Vpo 

Substituting this in (A4) gives 

+1 (^Biu) + l{h')B"{u) - {h')MB'{u)^ V {A8) 

Note that this equation has the same general form as 
the original equation (Al). If we denote the matrices 
A and B in the original equation (Al) as A{uX) and 
B{u, C) to indicate that these where the matrices ob- 
tained after eliminating all wavevector components of h 



with q > Cgo, then the new matrices obtained by ehm- 
inating the additional roughness with wavevectors be- 
tween {( — A^)g'o < q < C<lo becomes 

A{u, C - AC) = A{u,0 + ^{h^) A" {u,C) 



{h')A'iu,OMA'{u,C) 



{A9) 



or 



qq 

Cdct>C{ci) 



where |q| = C^o- Thus to first order in (/i^)^: 
A{u, 0=u^ + u3(l - 3D) 



(Alb) 



B{u, C - AC) = B{u, + - {h-')B"{u, Q 



-{h^)A'(u,C)MB'{uX) 



w^l 1 



3(1 -3£)) 



(A16) 



(AlO) 



Since B ^ u and — >■ as C — ^ Ci we can write 



Since is smah we can expand the left hand side to 
linear order in Ai^. Furthermore note that 



AC 



AcZ 



(C-AC)go<g<C9o 



d'q C(q) 



B{u, 0=u + 6i(u)(/i')f + b2{u){h^f^ + .. 
Substituting this in (A13) gives 

bi = -SMu^. 
Thus to first order in (/i^)^ we get 



{A17) 



f^"" dqq r dcf> C{qcos<f), qsincf)) B{u, 0=u- {h'')^u = u (l - ^^^d) 

ACj(C-AC)go ■''o \ " / 



^^^^3D] (A18) 



Jo 



C(Cgocos0, Cgosinf/)) 



= -^ f d'q dd) = ~{h\ {All) 

where is the ensemble averaged of the square of 

the roughness amplitude including only roughness with 
wavevector |q| > Cg'o- Thus from (A9), (AlO) and (All) 
we get 



dA 

ac 



A"{u,C)-A'iu,C)MA'iu,C) 



dC 



dB 
5C 



]^B"{u,Q-A'{u,QMB'{uX) 



{h\ {A12) 



(AIS) 

If we assume that D{() is independent of C, it is easy to 
solve these equations using perturbation theory to arbi- 
trary order in the surface roughness amplitude h. Since 
A^ and — >■ as C — >■ Ci we can write 



A{u, C)=u^ + ai(w)(/i^)f + a2{u){h'^)l + .. 
To first order in we get from (A12) 

ai = 3m - 9Mu^ 

where 



{AU) 



M = 



It is strait forward to calculate the higher order terms 
(e.g., a2 and 62) in the expansions (A14) and (A17) but 
here we will only do so for the case of surface roughness 
with isotropic statistical properties. In this case Aij = 
A{u,()Sij and Bij = B{u,QSij- Thus the matrix M in 
(A7) becomes 

Mij = A-\hY^ J d^q C{<i)q-^qiqj = ^^"1% 



and (A12) and (A13) reduces to 



dA _ 1 
5C ~ 2 



A"{u,0 



[A'{u,or 



A{u,0 



d_ 
dC 



{A19) 



dB 
5C 



B"{u,C) 



A'{u,C)B'{u,C) 



A{u,C) 



dC 



where A and B arc now scalar fields. Substituting (A14) 
in (A19) gives to second order in (/i^)^ 

ai + 2o2(/i^)c = + K - 6w-^a'i + Qaiu"^) 

or 



ai = --U 



1 9 
02 = - (a'/ — 6u~^a\ + 9aiu~^) = — 
4 ^ -"8 



Thus, to second order 



A=u'-Uh\-iu-\hri 



-u^d 

~ y 2 8 J 

In a similar way one obtain to second order 
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